First Principles Study of Zn-Sb Thermoelectrics 
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We report first principles LDA calculations of the electronic structure and thermoelectric prop- 
erties of /3-Zn4Sb3. The material is found to be a low carrier density metal with a complex Fermi 
surface topology and non-trivial dependence of Hall concentration on band filling. The band struc- 
ture is rather covalent, consistent with experimental observations of good carrier mobility. Calcu- 
lations of the variation with band filling are used to extract the doping level (band filling) from the 
experimental Hall number. At this band filling, which actually corresponds to 0.1 electrons per 22 
atom unit cell, the calculated thermopower and its temperature dependence are in good agreement 
with experiment. The high Seebeck coefficient in a metallic material is remarkable, and arises in 
part from the strong energy dependence of the Fermiology near the experimental band filling. Im- 
proved thermoelectric performance is predicted for lower doping levels which corresponds to higher 
Zn concentrations. 
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There has been a recent revival of activity in the search 
for improved thermoelectric |inaterials, with an empha- 
sis on new materials systemstl. The efficiency of a ther- 
moelectric is characterized by a dimensionless figure of 
merit, ZT = aS'^T/n = S'^/L, where T is the temper- 
ature, a is the electrical conductivity, S is the Seebeck 
coefficient (thermopower) and k is the thermal conduc- 
tivity, which contains both electronic and lattice contri- 
butions, K — Kci -\- Kiat- The ratio L = k/uT, which 
is often called the Lorentz number, is ordinarily lim- 
ited from below by its electronic value, Koi/ctT, given 
by the Wiedemann-Franz value, L — (7r^/3)(fcB/e)^. 
Current thermoelectric materials have ZT « 1. With 
the Wiedemann-Franz value for L, ZT > 1 requires 
S > 160/xV/K. Since doped semiconductors naturally 
have large thermopower, much of the thermoelectric re- 
search over the past 40 years has focused on covalent 
semiconducting compounds and alloys, composed of 4th 
and 5th row elements, with a view to finding low ther- 
mal conductivity materials that have reasonable carrier 
mobilities and high band masses, e.g. Bi2Te3, Si-Ge and 
PbTe compounds. Despite the research efforts spanning 
three decades, little progress in increasing ZT has been 
achieved until recently, and in particular Bi2Te3/Sb2Te3 
has remained as the material of choice for room temper- 
ature applications. 

In the last two years, howeverp. three new materials 
with ZT > 1 have been reportedoci, and these do not 
clearly fall into the same class as previous thermoelectric 
materials. One of these novel materials, (i-Zn^Shs, with 



reported ZT « 1.3 at temperature relevant to waste heat 
recovery, has a large region of linear temperature depen- 
dence of the resistivity (Fig. |l|) , suggestive of a metallic 
rather than semiconducting material. But unlike normal 
metals this is accompanied by high thermopowers. 
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FIG. 1. Experimental resistivity of /3-Zn4Sb3. Replotted 
with a linear temperature scale using data of Ref . [2] . 
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In this report, we present first principles calculations, 
within the local density approximation (LDA), simi- 
lar to our pjceijious calculations for binary and filled 
skutteruditesETu. Our main purpose is to understand the 
rather remarkable thermoelectric properties of /3-Zn4Sb3 
by analysing its transport properties based on electronic 
band structures. 




FIG. 2. Crystal structure of /3-Zn4Sb3. White and dark 
grey spheres denotes the 12 Zn and 4 Sb atoms on pure sites 
in the 22 atom rhombohedral unit cell. 6 light grey spheres 
are occupied by 89% of Sb and 11% of Zn atoms. 



site reported to have approximately 11% of Zn and 89% 
of Sb (light grey spheres in Fig. ||) in the actual structure 
was taken to be a pure Sb site. This adjustment yields 
a formula Zn^Sbs with 22 atoms per rhombohedral unit 
cell. The mixed occupancy was accounted for in a rigid 
band model as discussed below. Of the five Sb atoms 
three occupy the mixed site (Sb^™)), and two reside in a 
pure Sb site (Sb*-^-* ) which is crystallographically inequiv- 
alent. Two Sb'-^-' atoms form Sh^^ dimers parallel to the 
rhombohedral axis. The bond length between these two 
Sb atoms is 2.82 A, which is exactly twice of Sb covalent 
radius. This suggest covalent Sb'-^-'-Sb^^-' bonding, which 
is confirmed by the band structure calculations. Each 
Sb(P) forms three additional bonds with Zn, thus being 
essentially four-fold coordinated. The partial densities of 
states (Fig. ||) shows that these bonds are also largely 
covalent. 
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The computations were performed self-consistently 
within the LDA using the general potential—linearized 
augmented planewaves (LAPW) method@'E£l. This 
method makes no shape approximations to either the 
potential or charge density and uses flexible basis sets 
in all regions of space. As such it is well suited to ma- 
terials with open crystal structures and low site sym- 
metries like /3-Zn4Sb3. Special care has been taken to 
obtain a well converged basis set of approximately 2100 
functions with LAPW sphere radii of 2.5 a.u. for both 
Zn and Sb atoms. During the self-consistency iterations, 
the Brillouin-zone integrations were carried out with 10 
special k points in the irreducible wedge of the Brillouin 
zone; 781 inequivalent k points were used to calculate 
Fermi surface averages. Use of this relatively large num- 
ber of k points was necessary to obtain accurate Fermi 
velocities, and, especially, reliable Hall coefficients. We 
used the local orbital extensiontil of LAPW method to 
relax linearization errors in general and to include the 
upper core states consistently with the valence states. 
The calculations were rbased on the experimental crystal 
structure of /3-Zn4Sb3ll3 shown in Fig. ||; however, the 



FIG. 3. Total and partial density of states of /3-Zn4Sb3. 
Total density of states is for 22 atom unit cell and partial den- 
sity of states are for each inequivalent atom and magnified by 
a factor 10 for clarity. Note the highly covalent character of 
the band structure. The calculations were performed for the 
stoichiometric composition {Ep at lower vertical line). The 
actual Fermi energy is given by the upper marker. 

Both the Zn and Sb^'") occur in high coordinated po- 
sitions; Zn forms one bond with Sb(p), one bond with 
another Zn, and 3 bonds with Sb^™^ while Sb*^™-* forms 
all six bonds with Zn. The similar coordination of Zn and 
Sb(™) favors for Zn to substitute on Sb^") over Sb^P) , as 
observed. Substituting one Sb'™-' by a Zn eliminates 6 
Zn-Sb bonds and creates 6 Zn-Zn bonds. In general, this 
is not an energetically favorable process, and one could 
expect stoichiometric compound ZngSbs rather than Zn- 
enhanccd Zn6.33Sb4.77, as observed experimentally. On 
the other hand, in the stoichiometric compound, as dis- 
cussed below, the bonding states of Sb are not yet com- 
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pletely filled, and doping electrons into the system pro- 
duces an additional bonding effect. In other words, sub- 
stituting Sb by Zn reduces the number of Zn-Sb bonds, 
but makes the remaining bonds stronger. The balance 
between the two effects may account for the partial sub- 
stitution, however, the very narrow range of compositions 
for which samples exist remains unexplained in terms of 
the calculations. 

As mentioned, our calculations were done for the stoi- 
chiometric structure. The Fermi level lies « 0.4 eV below 
the top of the valence band, which is separated by a size- 
able gap (also « 0.4 eV) from the conduction bands. The 
effect of the deviation from stoichiometry can be taken 
into account in a rigid band model: it is assumed that the 
electronic structure of the actual compound is the same 
as for the stoichiometric one, the only difference being 
position of the Fermi level. This is expected to be valid 
in a highly covalent, broad band system as we find, and 
is supported a posteriori by the comparison with exper- 
iment. Correspondingly we have calculated the relevant 
transport properties as functions of the Fermi level posi- 
tion, and used the experimental Hall conductivity to find 
the Fermi level position corresponding to actual samples. 

The relevant transport properties were determined 
from the calculated band structures using the-standard 
kinetic theory as given by Ziman and otherdijO. For 
the electrical conductivity, the Bloch-Boltzmann kinetic 
equation in lowest order along a;-direction is (similarly 
for the other Cartesian directions) 



aAT)^e' / deN{e)vl{e)Tie,T) 



(1) 



Here N{e) is the density of electronic states at the energy 
e per unit volume, r is the scattering rate for electrons, 
and the quantity v^(e) is defined by: 
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The integrations are carried out over the iso-energy sur- 
face, dS'e, defined by ek = e. For = Ep, the Fermi 
energy, the integral is over the Fermi surface and is 
the Fermi velocity. In this case it is related to the square 
of the plasma frequency lo^^ = ATre'^N{Ep)v'^{EF), 
which is proportional to the optical carrier concentration 
(n/m)off — N{Ep)v-^{Ep). We note that in general this 
does not have a simple relationship to the Hall concen- 
tration or the electron count (i.e. the doping level); the 
Hall concentration is a measure of the average curvature 
of the Fermi surface, while the doping level is determined 
by the volume enclosed by the Fermi surface. 

For isotropic scattering, which is often a reasonable 
approximation, the relaxation time does not. cjrLtpr the 
expression for the Hall concentration, yieldingJlSO: 



(4) 



with 
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where for simplicity we have given the expression for 
cubic symmetry. Here M is the k-dependent effective 
mass tensor, defined as 



(6) 



Provided that the Fermiology does not vary strongly 
on the scale of k-gT, the derivatives of the Fermi distri- 
bution in the above expressions may be replaced by the 
T = limit, which is the delta function, thereby sup- 
pressing the energy integrals. This approximation was 
used in the calculations of the Hall concentrations in this 
paper. However, in the expression for the Seebeck co- 
efficient, 5", one has to include the energy dependence 
explicitly at temperatures several times smaller than the 
characteristic electronic energy scale. Thus we have used 
the full expression, 



Sm^^^jdee^T) 



dm 



de 



(7) 



where 



a{e,T)^e'N{e)vl{e)mT) 



(8) 

is the conductivity corresponding to a Fermi level posi- 
tioned at e. The expression for a{e,T) in Eq. (H) con- 
tains two energy dependent factors; a factor related to 
the square of plasma frequency, = Aire'^Nv'^ and the 
relaxation time, r. Ordinarily, the first ternipis the most 
energy dependent, but there are exceptionalj, e.g. Pd 
metal where Ep occurs near a very sharp feature in NLe^ 
and Kondo systems where there is resonant scatteringlZI. 
Since /3-Zn4Sb3 does not show any indication of such be- 
havior, we have approximated the energy dependence of 
a using only the Nv'^ term. 

P-Zn^Sh^ has a small Hall number and is character- 
ized either as a low carrier density metal or a heav- 
ily doped semiconductor. The former terminology is 
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probably more appropriate, because the resistivitya in- 
creases linearly with T over a wide range of at least 
300 K < T <550 K (See Fig. 0). The Hall number, 
uh = 9 X lO^^cm^"^ — 0.05 holes/cell, for the reported 
high ZT sample, although small for a metal may be too 
large to allow analysis of the transport in terms of usual 
semiconductor formulae. The Hall concentration calcu- 
lated according to Eq. (^) is plotted in Fig.^ as a function 
of the Fermi level shift from its position in the stoichio- 
metric compound. Note that uh is not simply related 
to the actual number of holes; the experimental Hall 
number of 0.05 holes/cell corresponds to the doping level 
of 0.1 holes/cell, roughly twice nn- Finally, the calcu- 
lated plasma frequency at the same position of the Fermi 
level is = 1 ~ 1.2 eV (depending on polarization), 
which corresponds to optical effective carrier concentra- 
tion {n/m)cs ~ 0.5(r7io/™) holes/cell with mp being the 
bare electron mass. Analysis of the Fermi surface (see 
Fig. shows that close to the top of the valence band 
in the range relevant to high ZT samples the main sheet 
of the Fermi surface is more toroidal rather than ellip- 
soidal. Thus even at this low carrier density, the Fermi 
surface has a complicated shape with electron-like and 
hole-like contributions to the Hall coefhcient. This ex- 
plains why the Hall concentration is so much different 
from the doping level, and also suggests that the Fermi 
velocity and conductivity may depend on the Fermi level 
position in a strong and unusual way. This is the case 
and it is partially responsible for the high thermopower. 
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FIG. 5. Main Fermi surface of /3-Zn4Sb3 at the experi- 
mental band filling. 
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FIG. 4. Calculated Hall concentration as a function of the 
Fermi level position relative to the stoichiometric solid (solid 
line) in a linear- log plot. The inset shows the same quan- 
tity over wider range including both sides of the gap. The 
dotted horizontal line represents the experimental value and 
the dot-dashed vertical line denotes the required Fermi level 
shift to match the calculated Hall concentration to measured 
value. 



FIG. 6. Calculated (short dashes) and experimental (solid 
line with open circle) thermopower for /3-Zn4Sb3 with 0.05 
hole/cell Hall concentration (0.10 hole/cell doping). S{T) 
for Hall concentrations 0.12 and 0.02 hole/cell (0.15 and 
0.04 hole/cell doping level) are given by the curves above 
(dot-dashed) and below (long dashed) of those for the actual 
doping level for comparison. 

Indeed, using Eq. (|^), we find rather high Seebeck co- 
efficient (Fig. H). Naturally, S depends strongly on the 
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position of the Fermi level. Comparison with the experi- 
ment shows that very good agreement is achieved for the 
the hole count of 0.1 hole/cell. This is exactly the hole 
count that we deduced from the measured Hall number 
of 0.05 holes/cell, demonstrating the consistency of our 
approach. 

As mentioned, the experimental resistivity up to ap- 
proximately 550 K shows typical metallic behavior with 
high-tcmpcrature resistivity coefficient dp/dT w 3.1 
/xfi-cm/K and residual resistivity po ~ 1.8 mf2-cm. From 
these data and the calculated plasma frequency one can 
estimate the transport clcctron-phonon coupling con- 
stant Atr and the scattering rate 7 due to static defects 
to obtain Atr ~ 1, and 7 » 0.2 eV. Both numbers are rel- 
atively large, but may represent overestimates since they 
are not based on single crystal data. In any case, the 
experimental data indicate that at high temperature k 
is mainly k^x in the best reported sample. This implies 
that variations in doping level that increase S will also in- 
crease ZT . Our calculations show not unexpectedly that 
raising the Fermi energy corresponding to lower carrier 
concentration, leads to higher values of 5, particularly at 
intermediate temperatures ranging from room tempera- 
ture to the maximum. This corresponds to higher Zn 
concentrations on the mixed site. 

In summary, we have studied the band structures and 
transport properties of /3-Zn4Sb3 using first principles 
electronic structure calculations. Acciording to our ab ini- 
tio calculations, the Fermi surface has complicated topol- 
ogy and this makes the estimation of doping level non- 
trivial. Combining our calculations with the measured 
Hall number for the actual material, we arrive at a rela- 
tively large, essentially metallic, value for carrier concen- 
tration, 0.1 hole/cell. Using this, we calculated Seebeck 
coefficient and found it to be in excellent agreement with 
the experiment both in absolute value and temperature 
dependence. Based upon these results, we characterize 
this material as metal with complex, energy-dependent 
Fermi surface, which provides large thermopower for rel- 
atively high carrier concentration. This may be useful 
in identifying other candidate thermoelectric materials. 
The /3-Zn4Sb3 system is not yet optimized for thermo- 
electric application and further increase in ZT is ex- 
pected if the doping level can be reduced. One important 
issue revolves around the question of why it is so difficult 
to make samples with varying Zn concentrations on the 
mixed site, since higher Zn concentration would lead to 
higher ZT. We speculate that the reason has to do with 
competition from other phases during the high temper- 
ature synthesis. It would be very interesting to attempt 



growth by more non-equilibrium processes such as pulsed 

laser deposition. 

The authors are grateful for helpful discussions with 
T. Caillat and J. P. Fleurial. Computations were per- 
formed using DoD HPCMO faciUties at NAVO and ASC. 
This work is supported by DARPA. Work at NRL is sup- 
ported by ONR. 



^ G. Mahan, B. Sales and J. Sharp, Physics Today 50 #3, 
42 (1997). 

^ T. Caillat, J. P. Fleurial and A. Borshchevsky, Proceed- 
ings ICT'96, edited by T. Caillat, A. Borshchevsky and 
J. P. Fleurial (IEEE Press, Piscataway, 1996), p. 151. 

^ B.C. Sales, D. Mandrus and R.K. Williams, Science 272, 
1325 (1996). 

" J.P. Fleurial, A. Borshchevsky, T. Caillat and G.P Meisner, 
Proceedings ICT'96, edited by T. Caillat, A. Borshchevsky 
and J.P. Fleurial (IEEE Press, Piscataway, 1996), p. 91. 

^ D.J. Singh and I.I. Mazin, Phys. Rev. B 56, R1650 (1994). 

®L. Nordstrom and D.J. Singh, Phys. Rev. B 53, 1103 
(1996). 

''D.J. Singh and W.E. Pickett, Phys. Rev. B 50, 11235 
(1994). 

® J.L. Feldman and D.J. Singh, Phys. Rev. B 53, 6273 (1996); 

54, 712 (1996) 
^ O.K. Anderson, Phys. Rev. B 12, 3060 (1975). 

D.J. Singh, Planewaves, Pseudopotential and the LAPW 

Method (Kluwcr, Boston, 1994). 
" D.J. Singh, Phys. Rev. B 43, 6388 (1991). 
^2 L. Hedin and B.I. Lundqvist, J. Phys. C 4, 2064 (1971). 
^'^ P. Villars and L.D. Calvert, Pearson's Handbook of Crys- 

tallographic Data for Interrnetallic Phases, 2nd ed. (ASM 

International, Materials Park, OH, 1991), p. 5206. 

G.A. Slack, in Solid State Physics, edited by F. Seitz and 

D. TurnbuU, Vol. 34 (Academic Press, New York, 1979), p. 

1. 

D.G. Cahill, S.K. Watson and R.O. Pohl, Phys. Rev. B 46, 
6131 (1992). 

J.M. Ziman, Principles of the Theory of Solids (Cambridge 
University Press, Cambridge, 1972). 
^"^ CM. Hurd, The Hall effect in Metallic and Alloys (Plenum, 
New York, 1972). 

W.W. Schulz, P.B. Allen and N. Trivedi, Phys. Rev. B 45, 
10886 (1992). 

A.E. Karakosov, I.I. Mazin and Y.A. Usperiski, Sov. Phys. 
Doklady 277, 848 (1984). 



5 



